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Waking the Colored Plasma 
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We calculate the wake induced in a hot, dense QCD medium by a fast parton in the framework 
of linear response theory. We discuss two different scenarios: a weakly coupled quark gluon plasma 
(pQGP) described by hard-thermal loop (HTL) perturbation theory and a strongly coupled QGP 
(sQGP), which had the properties of a quantum liquid. 


I. INTRODUCTION 

The quenching of QCD jets created in relativistic 
heavy ion collisions has been proposed as an important 
indicator for the creation of a quark-gluon plasma dll 
It is extensively studied theoretically and experimentally 
at RHIC. 

The main emphasis in the theoretical studies has been 
on the description of the radiative energy loss which 
the leading parton suffers when traveling through the 
medium due to the emission of a secondary partonic 
shower. For an overview of the quantum held theoret¬ 
ical description of this process and its theoretical and 
experimental implications, see e. g. Hill 

In this work we focus on another aspect of the in¬ 
medium jet physics, namely the fact that a charged par¬ 
ticle with high momentum traveling through a medium 
induces a wake of current density, charge density and 
(chromo-)electric and magnetic held conhgurations. This 
wake rehects signihcant properties of the medium’s re¬ 
sponse to the jet. 

We apply methods of linear response theory to the sys¬ 
tem of a relativistic color charge traveling through a QCD 
plasma. To investigate the coherence behavior of the 
plasma reacting to a charged particle traveling through 
the plasma, we consider the plasma’s reponse to that ex¬ 
ternal current. We restrict ourselves to only one weak, 
external current which points in a hxed direction in color 
group space, so that the held strength is linear in the cur¬ 
rent. In this framework quantum and non-abelian effects 
are included indirectly via the dielectric functions, and 
ct- Linear response theory implies that the deviations 
from equilibrium are small enough that they will not, in 
turn, modify the dielectric function. Furthermore, for 
simplicity, we assume a homogenous and isotropic exten¬ 
sion of the medium and disregard any hnite size effects. 

We discuss two qualitatively different scenarios: In the 
hrst one, we assume the plasma to be in the high tem¬ 
perature regime, where the gluon self-energies can be de¬ 
scribed using the leading order of the high temperature 
expansion, T 3> w, fc. In the second scenario we investi¬ 
gate what happens if the plasma behaves as a strongly 
coupled quark gluon plasma (sQGP) which can be de¬ 


scribed as a quantum liquid. Our consideration of the 
second scenario has been motivated by earlier studies of 
wakes induced by fast electrically charged projectiles in 
the electron plasma of metallic targets [y, 13 j well as 
by recent work exploring the induction of a conical flow 
by a jet in a sQGP Q. This latter study invokes a hy- 
drodynamical description of the energy deposited by a 
quenched jet in the medium and emphasizes the emer¬ 
gence of a Mach cone. The idea of determining the sound 
velocity of the expanding plasma from the emission pat¬ 
tern of the plasma particles traveling at an angle with 
respect to the jet axis has also been mentioned in j^. 

II. PLASMA PHYSICS 

The formalism for the linear response of the plasma 
can be straightforwardly generalized from electromag¬ 
netic plasma physics (see e. g. M)- In this formulation, 
a dielectric medium is characterized by the components 
of the dielectric tensor eij{uj,k), which for an isotropic 
and homogenous medium can be decomposed into the 
longitudinal and transverse dielectric functions cl and 
ct via: 

Gj = (1) 

Here = kikj/k'^ and Vt = 1 — Vl are the longitudi¬ 
nal and transverse orthonormal projectors with respect 
to the momentum vector k. One can relate the dielectric 
functions e l and €t to the self-energies H^ and Hy of the 
in-medium gluon via m-- 


1 1,\ 1 ^L{uj,k) 

ei(w,fc) = l ^2 ’ 

(2) 

t 1 nT(w,fc) 

eT[uj,k) = l 

( 3 ) 


where = uj'^ — k"^. 

Using Maxwell’s equation and the continuity equation 
in momentum space the total chromoelectric field in 
the QCD plasma is related to the external current jext 
via: 

ei’Pi + Vt 


Etoti^,k) = -r^j^^^{uj,k). (4) 
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ishes: 


det 


cl'Pl + (ct 


j)'Pt 


= 0 . 


(5) 


This equation governs the dispersion relation for the 
waves in the medium. Since it can be diagonalized into 
purely longitudinal and transverse parts, dispersion rela¬ 
tions for the longitudinal and transverse dieelectric func¬ 
tions can be derived m-- 

CL = 0, (6a) 

ex = {k/uj)^. (6b) 


These equations determine the longitudinal and trans¬ 
verse plasma modes. The longitudinal equation is also 
the dispersion relation for density-fluctuations in the 
plasma, namely space-charge fields which could be spon¬ 
taneously excited in the plasma without an application 
of external disturbances [1(T| . 

The color charge density induced in the wake by the 
external charge distribution is: 


Pind = 


{ex 


(7) 


The induced color charge density can also be calculated 
from the induced scalar potential via a Poisson equation: 


47r 

^ind — Pind 5 


( 8 ) 


Casimir invariant C (which is either Cf = 4/3 for quarks 
or Ca = 3 for a gluon). In this simplified model one 
disregards changes of the color charge while the particle 
is propagating through the medium by fixing the charge’s 
orientation in color space [Hill 

The non-radiative part of the energy loss of the inci¬ 
dent color charge is given by the back reaction of the 
induced chromo-electric field onto the incident parton. 
The energy loss per unit length is given by [l^. 111 : 

^ (13) 


where the induced chromo-electric field is the to¬ 
tal chromo-electic field minus the vacuum contribution. 
Using the inverse of Eqn. 0 it is given by: 


rpa _ 

-^ind — 


— - l) Vl + 
ei / 

1 


1 


ex - -xi 


1 - 


fc2 


Vx 


47r - 


-Jexf 


(14) 


From and la, the non-radiative energy loss per unit 
length is given by H: 


dE Cas 

dx 2 'k'^v J 

(z;2fc2 -u;2)Im(w2 


d^k^{lme 


L 

2\-l\ 




(15) 


if one works in a gauge where the vector potential is 
transverse to the momentum [Hill 

Since one can relate the total chromo-electric field to 
the induced charge in linear response theory by 

J/nd = *^^(1 - e)4“ot/(47’-), (9) 

a direct relation between the external and the induced 
current can be derived using Eqn. 0 : 


^ind 


— - 1 ) Ul + 

CL 


1 — Ct 


er — 


A;2 


'Ft 


Jext ’ 


( 10 ) 


The induced charge and the induced current obey a con¬ 
tinuity equation: 


'^k • Jind 7Cdpind — 0. 


( 11 ) 


where uj = v ■ k. 


III. CHARGE WAKE IN A QGP IN THE HIGH 
TEMPERATURE APPROXIMATION 

In this section we study the first scenario in detail. The 
plasma is assumed to be in the high temperature regime, 
where the gluon self-energies can be described by the so- 
called hard-thermal loop (HTL) approximation, the lead¬ 
ing order of the high temperature expansion T ^ uj,k 
[Hlill- Tlii® regime can be expected to be realized above 
the deconfinement temperature Tc- The self-energies de¬ 
rived in this approximation have been shown to be gauge 
invariant m and the dielectric functions obtained from 
these are therefore also gauge invariant. The dielectric 
functions read explicitly: 


At this point we specify the current and charge densi¬ 
ties associated with a color charge as Fourier transform 
of a point charge moving along a straight-line trajectory 
with constant velocity v: 

J'ext = 2.TTq°'v5{uj-v-k), (12a) 

P^xt = 2TTq°‘5{w-v-k) (12b) 



ITT© (l 



(16) 


where its color charge defined by = Cag with the where x — to/k. We first discuss the induced charge 
strong coupling constant as = and the quadratic and current densities pind and jind using eqs. (1711011 . The 
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Fourier transform of o reads in cylindrical coordinates: 


Pv,indip,z,t) = f / dn' k'Jo{k' pnig) X 

[Ziryv Jq 

J duj' exp iu)' ^ - rUg ^-, (17) 


ing that the induced charge density pv,ind is proportional 
to The cylindrical symmetry around the jet axis 

restricts the form of the current density vector j-md- It 
has only non-vanishing components parallel to the beam 
axis, jv,ind, and radially perpendicular to it, jp^md- 


where k = sjxP' ^ jv'^, k = n'mg and w = to'mg, show- 


J 


^3 poo pOO 

jv,ind{p,z,t) = I ^ q°- / dn' k'Jq{k' pmg) / duj'exp 

(Ztt) V Jq J_^ 


im 


jp,znd{p,z,t) = i q°- / dn' k!Ji{k' pmg) / dw'exp 

(27r)2u 7o J-oo 


Again, the components of the current density are propor¬ 
tional to m^g. These expressions are only valid for ?; > 0. 
For a parton at rest the current density vector vanishes 
and the induced charge density can be directly derived 
from Eq. 0 . It has a Yukawa-like shape and reads (for 
r = |71): 


Pv=o.ind(?') = --^exp{-\/2mgr)q°^. (19) 

ZTrr 

Since longitudinal and transverse plasma modes, as de¬ 
termined by the dispersion relations dini, can only ap¬ 
pear in the time-like sector of the uj,k plane [ij, Il3ll - 
collective excitations do not contribute to the charge and 
current density profile of the wake. Emission analogous 
to Cherenkov radiation and Mach cones do not appear, 
but the charge carries a screening color cloud along with 
it. Eig. n illustrates this physically intuitive result nu¬ 
merically. It shows the charge density of a colored parton 
traveling with v = 0.99c in cylindrical coordinates. 

Despite the fact that Mach cones do not appear in the 
charge density wake, the particle still suffers energy loss 
due to elastic collisions in the medium, which can be 
described by eq. This effect of energy dissipation 

has been studied in [TJ and we refer the reader to this 
work for details. The integrand in 03 contributes to 
the integral in the space-like region only, where \x\ < 1, 
and therefore does not get contributions from frequencies 
where collective plasma modes exist. This is consistent 
with the fact that such modes are not excited in the in¬ 
duced charge and current densities. 


IV. CHARGE WAKE INDUCED IN A 
STRONGLY COUPLED QGP 

In this section we consider another scenario and as¬ 
sume that the plasma is in a strongly coupled regime 



where it exhibits properties characteristic of a quantum 
liquid. We investigate the similarities and differences in 
the wake of a colored parton in such a quantum liquid 
with those obtained for the perturbative, ultrarelativis- 
tic plasma discussed in the last section. Since up to date 
there is a lack of theoretical methods for first principle 
calculations of the color response functions in a strong 
coupled QGP regime, we restrict our discussion in this 
section to a simple model. Nevertheless, we construct 
this simplified model in such a way that important gen¬ 
eral conclusions about the differences and similarities of 
a quantum liquid scenario and the scenario of the last 
section can be drawn. 

The most prominent difference - in accordance with a 
sQGP quantum liquid scenario - is the possibility of a 
plasmon mode that extends into the space-like region of 
the uj, k plane above some k. 

The sQGP paradigm suggests very low dissipation at 
small k (’’hydro modes”), but large dissipation at high 
k. Our assumption is that a critical momentum kc sep¬ 
arates the regimes of collective and single particle exci¬ 
tation modes in the quantum liquid where the dominant 
modes below kc are plasmon excitations and where dis¬ 
sipative single-particle response is negligible. Since we 
are predominantly interested in collective effects in the 
plasma here, we restrict our study to the region below 
kc- 

To be specific, we assume that the dielectric function 
of the strongly coupled plasma in the k < kc regime leads 
to a longitudinal dispersion relation of the form: 


WL = -|-a;2 , (20) 

where cUp denotes the plasma frequency and u the speed of 
plasmon propagation which is assumed to be constant. In 
accordance with (EUll we assume the following dielectric 
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FIG. 1: Spatial distribution of the induced charge density 
from a jet with fixed color charge in a high temperature 
plasma where the HTL approximation is applicable. The 
lower plot shows equi-charge lines in the density distribution. 
The density proHle is a cloud traveling with the color charged 
jet. 


function 


Cl — 1 + 




+ uj'i/2 


{k < kc) ■ 


( 21 ) 


Note that this differs from the classical, hydrodynamical 
dielectric function of Bloch ca , since the latter one is 
singular at small k and uj due to phonon contributions 
which cannot mix with colored plasmons. The dielec¬ 
tric function m is constructed in such a way that it 
allows to study one possible aspect of a quantum liquid 
scenario, namely a plasmon mode which extends into the 
space-like region of the w, k plane. The qualitative in¬ 
duced wake structure in such a quantum liquid scenario 
is general, namely a Mach shock wave structure for a su¬ 
personically traveling color source. In that respect the 
principlal findings of this and the next section can be 
expected to hold generally for a quantum liquid with a 
plasmon branch similiar to (ED independent of the exact 
form of the dielectric function. 

We assume a speed of plasmon propagation of u/c = 


l/-\/3. This differs from the speed of plasmon propaga¬ 
tion obtained in the HTL approximation in the small k 
limit M, namely ujc = -^3/5, and of a hadronic res¬ 
onance gas ujcK. \/(L2 |la[l3' We disregard possible 
damping effects for k < k^- This assumption implies a 
negligible imaginary part in the dielectric function of the 
strongly coupled QGP, i. e. no collisional energy loss in 
the collective regime. 

Determining the plasmon mode via Eq. ED reveals 
that in this scenario the mode is in the space-like region 
of the w, k plane for k > ojp/[s/l — v?) . Recall that this 
is different for the high-temperature plasma, where longi¬ 
tudinal and transverse plasma modes only appear in the 
time-like region, |x| = \oj/k\ > 1. In the quantum liquid 
scenario one can expect that the modes with low phase 
velocity |x| < u suffer severe Landau damping because 
they accelerate the slower moving charges and decelerate 
those moving faster than the wave [ifl [13 • A charge 
moving with a velocity that is lower than the speed of 
plasmon propagation can only excite these modes and not 
the modes with intermediate phase velocities u < |a;| < 1, 
which are undamped The qualitative prop¬ 

erties of the color wake are in this case analogous to 
those of the high temperature plasma, namely that the 
charge carries only a screening color cloud with it and 
Cherenkov emission and Mach cones do not appear. Us¬ 
ing eq. ED and restricting the integration area to the 
region k < k^ = 2ujp,^ one can illustrate this for a col¬ 
ored particle traveling with v = 0.55c < u. Fig. [21 shows 
the charge density cloud traveling with the colored par- 
ton. 

If the colored parton travels with a velocity v > u 
that is higher than the speed of sound in the medium, 
modes with an intermediate phase velocity u < |a;| < 1 
can be excited. The emission of these plasma oscillations 
induced by supersonically traveling particles is analogous 
to Cherenkov radiation. This leads to the emergence of 
Mach cones in the induced charge density cloud with the 
opening angle 


A® = arccos (u/u) (22) 

The effect is well known in solid state physics M and 
is analogous to the Mach cones induced by fast heavy ions 
in electron plasmas [mill The physics of shock waves 
in relativistic heavy-ion collisions has also been discussed 

in miiiiii. 

Figures|Di,b show the induced charge density for a col¬ 
ored particle traveling at u = 0.99c > u. Figs.^,b show 
the corresponding current densities. Here eqs. JTTim 
have been used, and the integration area has been re¬ 
stricted to the region k < kc = 2u!p. Numerical consis¬ 
tency has been checked also via the continuity equation 

ED- 


^ Note that the expressions in HQ do not correspond to k < kc, 
but K. < kc. 
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FIG. 2: Spatial distribution of a induced charge density from 
a jet with high momentum and a fixed color charge q°' that 
is traveling with v = 0.55c < u. The lower plot shows equi- 
charge lines in the density distribution. The density profile is 
that of a cloud traveling with the color charged jet. 


We emphasize that the existence of Mach cones is ex¬ 
pected in a plasma in general if the particle is moving 
faster than the speed of sound in the plasma and if the 
dispersion relation of the collective mode extends into the 
space-like region. The wake induced by a colored jet in 
such a setting leads to regions of enhanced and depleted 
charge density in the wake, which have the shape of Mach 
waves trailing the projectile. 


V. OBSERVABLE CONSEQUENCES 

It can be expected that the phenomenon of Mach cones 
should eventually lead to a directed emission of secondary 
partons from the plasma. This has its anology in solid 
state physics where Mach cones induced by fast heavy 
ions in an electron plasma lead to an emission spectrum 
of electrons carried within the wake HQ which have 
been studied experimentally [^ . 

If that scenario of a strongly coupled QCD plasma is 
realized in the matter created in an ultrarelativistic heavy 
ion collision, one could expect to observe these cones in 
the actual distribution of secondary particles associated 
with jets at RHIC, a signature also proposed in [^l^. 

Indeed preliminary data from the STAR experiment 
(see Fig. 1 in j^) (also PHENIX data can be expected 
to be available regarding this issue) show such effects 
in the background distribution of secondary particles in 



FIG. 3: (a) Spatial distribution of the induced charge density 
from a jet with high momentum and f ixed color charge that 
is traveling with v = 0.99c > u = ^/TJ3c. (b) Plot showing 
equi-charge lines in the density distribution for the situation 
in (a). 


the azimuthal angle Ac/). The peak near zero degree 
corresponds to secondaries from an outmoving jet, par¬ 
ticles from the companion jet result in a distribution 
with a clear maximum near Ac/) = tt for pp collisions, 
where no medium effects are present. In Au-Au colli¬ 
sions there is a distribution with two maxima around the 
Ac/) = TT position. These are located at A^ « tt ± I.l. 
One can argue 0 that such an effect could possibly be 
traced back to a Mach shock front traveling with the side- 
away jet leading to maxima in the distribution at about 
A$ = 7r±arccos(u/u). Given the confirmation of this ef¬ 
fect in the data, this would clearly indicate that the first 
scenario, viz. a plasma described by perturbative color 
dielectric functions, is not realized in the RHIC experi¬ 
ments. Furthermore, constraints of the properties of the 
plasmons in general could be deduced. 

The preliminary data indicate that a pronounced Mach 
cone might be reflected only in the secondaries of the 
quenched ]et. This can probably be attributed to the cir¬ 
cumstance, that the unquenched jets are emitted points 
too close to the plasma surface to generate a fully formed 
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FIG. 4: (a) Spatial distribution of the equi-value lines of the 
current densitie’s c omp onent parallel to a jet travelling with 
V = 0.99c > u = i/l/Sc in a sQGP. (b) Spatial distribution 
of the equi-value lines of the current densitie’s comp onen t 
perpendicular to a jet travelling with v = 0.99c > u = ^JTJic 
in a sQGP. 

wake behind them. 

The speed of plasma propagation could be also deter¬ 
mined. In fact if the maxima at A(/) « tt ± 1.1 are exper¬ 
imentally confirmed, it would correspond to u/c « VO.2 
0. It is interesting to note, that a study of the an¬ 
gular structure of the collisional energy losses of a hard 
jet in the medium would also support such an observa¬ 
tion: the incident hard jet’s scattering angle vanishes in 


the relativistic limit - leading to the jet’s propagation 
along a straight line - whereas the expectation value of 
the scattering angle 0 of a struck ’’thermal” particle is 
(0) « 1.04 which is close to 1.1. 

VI. CONCLUSIONS 

In this letter we discussed the properties of the charge 
density wake of a colored hard partonic jet traveling 
through a QGP plasma in the framework of linear re¬ 
sponse theory. We have studied two different scenarios, 
namely high temperature QGP at T ^ Tc described in 
the HTL approximation, and a description of a strongly 
coupled QGP (sQGP) behaving as a quantum liquid. We 
found that the structure of the wake corresponds to a 
screening color cloud traveling with the particle in the 
case of the high temperature plasma and in the case of a 
quantum fluid, if the velocity of sound in the plasma is 
not exceeded by the jet in the latter case. The structure 
of the wake is changed considerably in comparison to the 
former cases, if the jet’s velocity exceeds the plasma’s 
speed of sound and the collective modes have a disper¬ 
sion relation extending in the space-like region. Then 
the induced parton density exhibits the characteristics of 
Mach waves trailing the jet at the Mach angle. 

It is argued that this effect could be used to constrain 
theoretically possible scenarios by experimental analy¬ 
sis via the angular distribution of the secondary parti¬ 
cle cones of jet events in RHIC. First indications of the 
observation of a Mach shock phenomenon in RHIG in 
the quenched jet’s secondary particle distribution from 
STAR data were discussed. 

In general secondary particle distributions can be used 
to provide methods of probing the QCD plasma’s collec¬ 
tive excitations experimentally. 
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